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Liu et al. [Phys. Rev. E 74, 056302 (2006)] have presented a WKB analysis of the helical magne- 
to-rotational instability (HMRI), and claim that it does not exist for Keplerian rotation profiles. We 
show that if radial boundary conditions are included, the HMRI can exist even for rotation profiles 
as flat as Keplerian, provided only that at least one of the boundaries is sufficiently conducting. 

PACS numbers: 47.20.-k, 47.65.-d, 52.30.Cv, 52.72,+v 



The helical magnetorotational instability (HMRI) [J, 
S H, 0, [H , is similar to the standard magnetorotational 
instability (SMRI) 0,0,1, in the sense that both are 
mechanisms whereby hydrodynamically stable differen- 
tial rotation profiles may be destabilized by the addition 
of magnetic fields. However, the way in which they be- 
have in the limit of small magnetic Prandtl number Pm 
is very different. Unlike the SMRI, which ceases to ex- 
ist for zero Pm, the HMRI continues to exist, with the 
relevant nondimensional parameters being the hydrody- 
namic Reynolds number Re and the Hartmann number 
Ha, but the magnetic Reynolds number Rm = Pm • Re 
and the Lundquist number S = Pm 1 / 2 • Ha (the relevant 
parameters for the SMRI) tending to zero along with Pm. 

Given these very different scalings, one natural ques- 
tion to ask is whether they exist for the same range of 
rotation profiles. The SMRI is known to operate for any 
outwardly decreasing profile^ In contrast, the HMRI is 
more delicate; the results of [l| indicate that as one moves 
beyond the Rayleigh line to ever flatter profiles, one even- 
tually switches back to the SMRI scalings. Liu et al. @ 
explored this issue more systematically, and claim that 
the HMRI cannot exist for profiles as flat as Keplerian. 
Specifically, their WKB analysis indicates that if Q. ~ r™ , 
then the HMRI requires n < —1.66, thereby excluding 
the Keplerian value n = —1.5. [Their Eq. (12), with 
their Ro = n/2.] 

We do not dispute the validity of their analysis; indeed, 
some of our results below are in excellent agreement with 
it. However, we note that any WKB analysis is neces- 
sarily local, and does not incorporate the boundary con- 
ditions of the global problem. We show here that if the 
radial boundary conditions are taken into account, the 
HMRI can exist even for the Keplerian value n = —1.5, 
provided only that one of the boundaries is at least some- 
what conducting. 

The eigenvalue problem to be solved is 

Rejv = D 2 v + Reikr- 1 (r 2 nyip + R& 2 ikb, 
Re 7 D 2 ifj = L>V - Re 2ik fl v + Ha 2 (k 2 ^ + 2ikf3r~ 2 b), 
= D 2 b + ikv- 2ik f3r~ 2 ip, 



essentially the same as in [IJ, except that we restrict at- 
tention to the Pm — > limit, ensuring that any instabil- 
ities we obtain will necessarily be the HMRI. (Priede et 
al. [l(| have very recently also considered this Pm — > 
limit; some of their results are quite relevant here, and 
will be discussed below.) 

The boundary conditions associated with v and tp are 
no-slip, just as in [l|. The boundary conditions associ- 
ated with b are 

b = e(rb)' at n = 1, 6 = at r„ = 2. 

The outer boundary is therefore insulating, whereas the 
nature of the inner boundary depends on e: e = is 
insulating, e = 00 is perfectly conducting. Intermediate 
values correspond to a boundary consistin g of a thin layer 
of relative conductance e (see for example [ill [l2| for this 
thin boundary layer approximation in other contexts). 

f2(r) is the rotation profile whose stability is to be in- 
vestigated. In p], Q we considered Taylor-Couette pro- 
files of the form = c\ + c^jr 2 . To facilitate comparison 
with Liu et at, here we will primarily consider profiles of 
the form = r n . As we will see though, the two choices 
yield almost identical behavior. 

Figure 1 shows contour plots of the critical Reynolds 
number for the onset of the HMRI, as a function of n and 
/?, and optimized over the Hartmann number Ha and the 
axial wavenumber k. Turning first to the e = plot on 
the left, we note how increasing from 1 to 5 facilitates 
the instability, that is, allows it to exist increasingly far 
beyond the Rayleigh line at n = —2. Beyond j3 w 5 
though the Re c curves become largely independent of j3. 
And crucially, even the Re c = 10 6 curve asymptotes just 
to the left of the Liu et al. line at n = —1.66. These 
e = results are therefore in excellent agreement with 
their prediction that the HMRI only exists to the left of 
this line. 

However, if we now turn to the e = 00 plot on the right, 
we note that Re c = 10 4 already extends beyond the Liu et 
al. line, and Re c = 10 5 extends beyond the Kepler line 
n = —1.5. Simply switching the inner boundary from 
insulating to conducting is sufficient to allow the HMRI 
to operate even for Keplerian rotation profiles. Similar 
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FIG. 1: The left panel is for e = 0, the right for e = oo. In 
each case the solid curves are contours of log Re c as a function 
of n and /3, optimized over Ha and k. The dashed curves are 
the location where Re c is optimized over (5 as well. The dotted 
curves show how the 10 3 contours are shifted slightly to the 
right if the r™ profile is replaced by a Taylor-Couette profile 
having the same £l /Qi ratio. The dash-dotted lines denote 
the n = —1.66 boundary, beyond which the HMRI does not 
exist in the Liu et al. analysis. Note how this agrees very well 
with our results for e — 0, but not at all for e — oo. 




FIG. 2: log Re c as a function of n, optimized over k, Ha and 
/3. Note how e = does indeed appear to reach a vertical 
asymptote at the Liu et al. value n = —1.66, but e = 0.5 is 
already sufficiently large to reach the Kepler value n = —1.5. 



results are obtained if instead it is the outer boundary 
that is switched from insulating to conducting. Why the 
electromagnetic boundary conditions should have such a 



dramatic effect is not clear, but it is certainly well known 
in many other contexts, e.g. [l3|, [Hj]> that they can play 
a crucial role. 

Note also that exactly the same phenomenon illus- 
trated in Fig. 1 is already implicit in Fig. 2 of [lOj, where 
the HMRI exists up to fj, « 0.32 for insulating bound- 
aries, but up to /j, ~ 0.45 for conducting boundaries, 
where /i = Q /Cli. Translating from [i to n via fi = 2™, 
their results become n s=s —1.64 for insulating boundaries, 
versus n sa —1.15 for conducting boundaries. That is, the 
Keplerian value n = —1.5 is accessible with conducting 
boundaries, but not with insulating ones. 

At first sight it would appear that their value of 
n w —1.64 for insulating boundaries is already in con- 
flict with the Liu et al. limit n < — 1.66. In fact, this 
slight discrepancy is due to the difference between the 
ft — ci + C2 /r 2 profile used by [l(| , and the f2 = r" pro- 
file considered by Liu et al. Even if ci and ci are chosen to 
match fl[ = 1 and Q Q = 2 n , a Taylor-Couette profile will 
be somewhat steeper near n, and correspondingly some- 
what shallower near r a . By concentrating near the inner 
boundary, the instability can then operate at somewhat 
larger f2 Q /f2i values than for an r n profile. The dotted 
lines in Fig. 1 quantify this effect, showing the Re c = 10 3 
curves if the r n profile is replaced by a Taylor-Couette 
profile, with ci and ci chosen as indicated above. We see 
that qualitatively the two profiles yield exactly the same 
behavior, but that the Taylor-Couette profile extends to 
slightly larger values of n. It was precisely to avoid this 
effect, and thereby allow a direct comparison with Liu et 
al., that we chose here to concentrate on the r n profiles. 

The dashed lines in Fig. 1 show the locations along 
which Re c is optimized not only over k and Ha, but over 
(3 as well. Fig 2 shows Re c as a function of n along 
these lines, now including not just e = and oo, but 
also e = 0.5 and 1. We see that e = 0.5 is already large 
enough to reach the Kepler value n — —1.5. That is, if 
the conductance of the inner boundary is only 1/2 that 
of the fluid region, the HMRI already exists for rotation 
profiles as flat as Keplerian. 
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